arXiv:hep-ph/9706285v2 22 Aug 1997 


SLAC-PUB-7528 

hcp-ph/9706285| 

June, 1997 


Complete 0(af) Results for e + e —■» (j,Z) —> Four Jets [] 

Lance Dixon and Adrian Signer 

Stanford Linear Accelerator Center 
Stanford University 
Stanford, CA 94309 


Abstract 

We present the next-to-leading order (0(aA)) perturbative QCD predictions for e + e“ 
annihilation into four jets. A previous calculation omitted the 0(a terms suppressed 
by one or more powers of 1 / N'f, where N c is the number of colors, and the ‘light-by-glue 
scattering’ contributions. We find that all such terms are uniformly small, constituting 
less than 10% of the correction. For the Durham clustering algorithm, the leading and 
next-to-leading logarithms in the limit of small jet resolution parameter y cu t can be 
resummed. We match the resummed results to our fixed-order calculation in order to 
improve the small y CVL t prediction. 
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1 Introduction 


Electron-positron annihilation into jets provides an arena for studying quantum chromody¬ 
namics (QCD) that is free of initial-state uncertainties such as parton distribution functions. 
At the large center-of-mass energies achieved by SLC, LEP, and now LEP2, e + e _ annihilation 
is also relatively free of nonperturbative final-state effects, i.e. hadronization corrections. On 
the other hand, perturbative QCD corrections to jet rates can be very large. For example, 
the three-jet rate at the Z° pole receives a 20-30% correction [|l| at order a^. These next-to- 
leading-order (NLO) corrections are of course critical for obtaining a precise experimental 
measurement of a s from the three-jet rate and related 0(a s ) observables [|2], H|. 

More recently, the NLO corrections to e + e _ production of four jets were computed, and 
a correction of roughly 100% was found [|4| for most jet algorithms (when the renormal¬ 
ization scale was set equal to the center-of-mass energy). This computation omitted terms 
suppressed by one or more powers of 1/A)?, where N c is the number of colors in a general 
SU(N C ) gauge theory (N c = 3 for QCD). It also neglected the 1 light-by-glue scattering’ 
contributions — interference terms where two different flavor quarks couple to the virtual 
photon or Z boson. In this article we present the complete 0(a^) results, using an im¬ 
proved version of the same numerical program, MENLO_PARC ||, which was employed 
for the leading-in-W computation. The crucial ingredients for the construction of the pro¬ 
gram are the tree-level amplitudes for five massless final state partons, e + e _ —> qqggg and 
e+e” -*■ qqq'q 'g [| 0, and especially the recently-computed one-loop virtual amplitudes for 
e + e - — *• qqq'q' || |Di] and e + e _ —»• qqgg |nj. We use the formulas given in refs. BB0. 

The NLO prediction of the four-jet fraction — an observable whose expansion begins at 
order a^ — makes it possible to measure cc^ 3 with the same formal level of precision (NLO) 
as has previously been reserved for 0(a s ) observables in e + e _ annihilation. However, the 
theoretical uncertainty in such a measurement will still be sizeable: Because the one-loop 
corrections are so large, the renormalization-scale dependence of the NLO four-jet result is 
still strong, and it is likely that uncalculated higher-order corrections are important. Also, 
a significant four-jet rate only appears at smaller values of the jet resolution parameter y cut , 
where there are large perturbative logarithms, although these can be partially resummed for 
the Durham algorithm |TT] . 

There are at least two other motivations for studying e + e _ annihilation to four jets: (1) 
These events are a background to e + e _ —> W + W~ —> 4 jets, particularly when the center-of- 
mass energy is not far above the ID-pair threshold, as is the case at LEP2. (2) Four-jet final 
states provide QCD tests to which three-jet events are insensitive |12|]. For example, the 


non-abclian three-gluon vertex appears at leading order in four-jet events; the same is true 
for the production of hypothetical, light, colored but electrically neutral particles, such as 
light gluinos |L3], [14|, [L5]. [L6[]. In both applications, distributions of the four jets with respect 
to energies and angles [|l!| are important. Such distributions can be computed at NLO using 
the same numerical program, and will be the subject of a separate publication [|D]]; here we 
briefly study the sensitivity of the total four-jet rate to additional light fermions. 

The remainder of the paper is organized as follows. In section 2 we describe the de¬ 
pendence of the four-jet rate on electroweak and color factors, and outline the structure of 
the numerical calculation. In section 3 we present the complete 0(a^) predictions for three 
different jet algorithms. We indicate the dependence of the predictions on the (unphysical) 
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renormalization scale /j. The Geneva algorithm |18| has a relatively mild /i dependence 
(small NLO correction) and a relatively strong dependence on the number of light quark 
flavors Nf ; we discuss the extent to which Nf can be determined from the Geneva four-jet 
rate alone. In section 4 we present results from matching the resummed Durham jet rate 
to the fixed-order 0(a^) results; the improved prediction agrees quite well with preliminary 


SLD data 191. Section 5 contains our conclusions. 


2 Structure of the Cross-Section and Computation 

For computational reasons as well as to study the effect of varying parameters, it is useful 
to decompose the leading-order (Born) and NLO contributions to the four-jet differential 
cross-section with respect to both their electroweak and QCD (color) structure. To simplify 
the electroweak decomposition we assume that the observable being calculated is insensi¬ 
tive to both (1) correlations between the final-state hadrons and the electron-positron beam 
direction, and (2) quark and gluon hclicities. We also assume the positrons are unpolar¬ 
ized and the electrons have a longitudinal polarization of P e (P e = +1 for a right-handed 
beam). QED initial state radiation and other electroweak corrections are neglected. Then 
the hclicity-summed four-jet (differential) cross-section at center-of-mass energy y/s may be 
written 
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where a is the hne structure constant, Q q is the charge of quark q in units of e, and the left- 
and right-handed couplings to the Z° are 


V £ r = 


vi. = 


—1 + 2 sin 2 9 


w 


sin 29w 

±1 - 2 Q q sin 2 9 


V R ~ 


2 sin 2 9 


w 


w 


sin 29 


w 


Vr = 


sin 2 9w 
2Q q sin 2 9 W 
sin 2 9w 


(3) 


where 9\y is the weak mixing angle; the two signs in vf, correspond to up (+) and down (—) 
type quarks. Equations (0) and © include both virtual photon and Z boson exchange (and 
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their interference); the ratio of Z and photon propagators is given by 


Vz{s) = 


s — Ml + iT z M z ’ 


(4) 


where Mz and Tz are the mass and width of the Z. 

Representatives of the classes of diagrams contributing to f (n> and are de¬ 

picted in Fig. |l] as amplitude interferences. Five-parton cuts of these graphs, shown as 
dashed lines, correspond to the real part of the NLO correction; four-parton cuts, shown as 
dotted lines, correspond to the virtual part. In contribution (I) a single fermion couples to 
both ( 7 , Z) vector bosons in the interference, via either a vector or axial vector coupling. 
(As shown in the figure, there may be a second or even a third fermion loop in the interior 
of the graph, corresponding to ‘QCD’ factors of Nf in the cross-section.) This contribution 
dominates the cross-section at 0(a 2 ) and as we will see, again at 0(a^). 

The remaining contributions, (II) and (HI), have different origins in the real and virtual 
parts of the calculation. In the real part they come from the qqq'q'g final state when the 
roles of q and q' are exchanged on the opposite side of the cut; in particular, a different quark 
pair couples to the ( 7 , Z) on each side of the interference. In the virtual part they can have 
the same kind of exchange origin in qqq'q' final states, but they can also arise from either 
qqgg or qqq'q' graphs where a quark loop couples directly to the photon or Z (for example 
the contribution A 6;3 (Ag x ) in ref. [|]). 

Contribution (II) represents 1 light-by-glue scattering’, whereby a different fermion line 
couples to each vector boson, via a vector coupling in each case. There is no such contribution 
at 0(a 2 ) if only charge-blind observables are considered |2(J, due to Furry’s theorem — the 
order a 2 amplitude interferences all contain fermion triangle subgraphs. Although the cross- 
section is nonvanishing at 0(al), we shall see that it is still extremely small, due partly 
to cancellations in the sum over quark couplings in /R0(s), and partly to approximate 
cancellations in the phase-space integrations that are related to the exact cancellations at 
order a 2 . 

Contribution (III), ‘Z-by-glue’ scattering, is similar to contribution (II) except that the 
quarks couple to the Z through the axial vector coupling. This contribution is nonzero at 
although small for the three- and four-jet rates, and it remains small at O(a^). 
In Eq. (H) we have already carried out the sum over the five light quark flavors, in which the 
massless weak isospin doublets ( u , d ) and (c, s ) cancelled, leaving only the (t , b ) contribution. 
The top quark contribution to (III) is purely virtual for y/s < 2 m t , but it does not decouple 
in the large m t limit |£TJ . We expand in the limit of large top quark mass, including all terms 
through 0(s/m 2 )] at this order the top quark does not appear in the vector contribution 

(ii) aug. 

Dividing the four-jet cross-section <j 4 _j e t by the total hadronic cross-section at O^qls), 
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Figure 1: Representative contributions of type (I), (II) and (III), as described in the text. The 
coupling of a quark to the ( 7 , Z) vector boson is denoted by x, with a 1 ( 75 ) for vector (axial vector) 
coupling. Dashed lines correspond to representative five-parton cuts; dotted lines to four-parton 
cuts. 
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Neglecting for the moment the renormalization-scale dependence of the calculated cross- 
section we write the expansion in a s as 
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Next we decompose the one-loop correction to a 4 with respect to N c and Nf 


A = I, II, III, 
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Correspondingly, we write the full 0(a^) correction to the four-jet rate as 


/~t _ zA a ) I r^(P) I ri 
C ^4 — ^A l C4 -f C4 


(c) + ci d) + c 4 


(e) + ci f) + c 4 


(I!) + C (III) 


(9) 


( 10 ) 


we calculated c {°' ,b '^; 


The subleading-color terms C[ d ' e ’^ come 


absorbing all prefactors into the definitions of the C{ x> . In ref. 
here we add the remaining terms in Eq. ([TOj) 
partly from non-planar interference graphs (not shown in Fig. |1|). They include identical- 
quark Pauli exchange contributions analogous to the E terms of ref. [|20|, as well as various 
subleading-color virtual subamplitudes P. [HI, and subleading terms in the real and virtual 


color sums. We find that all the additional terms are considerably smaller than c[ a,b,c \ at 
least for the overall four-iet rate. 

(x) 

The Monte Carlo integrations required to numerically evaluate the C\ 2 are done sepa¬ 
rately for each term, except that and are combined. An advantage Q of breaking 
up the problem in this way is that the 1 / A^-suppressed integrands have significantly more 
complicated analytic representations than the leading terms, and therefore take more time 
per point to evaluate (in some cases up to a factor of five longer). On the other hand, the 
1/Aj? parametric suppression implies that far fewer numerical evaluations of the sub leading 
terms are required in order to achieve an absolute statistical accuracy comparable to that for 
the leading-in-7V C terms. Contributions (II) and (III) could have been further decomposed by 
analogy to Eq. (|D|), but in view of their small overall contribution they were each integrated 
as a single expression. 

As in any NLO QCD computation, the real and virtual corrections to the cross-section 
are separately divergent, but have a finite sum. I 11 dimensional regularization with D = 
4 — 2e, the singularities of the virtual part manifest themselves as poles in e in the one-loop 
amplitudes, whereas the real singularities are obtained upon phase-space integration of the 
squared tree amplitudes. We use a general version of the subtraction method |2(J to extract 
the singular parts of the real cross-section and combine them with the virtual poles. This 
method leaves a finite integral over five-parton phase space, and another over four-parton 
phase space, which are performed by adaptive Monte Carlo integration using VEGAS . 
The particular form of the subtraction method used here is essentially that described in 
ref. |2B[, to which we refer the reader for more details. No approximation of the matrix 
elements or the phase-space has to be made in this method. 


5 













The subtraction method relies on the fact that the integral over the tree cross-section is 
rendered finite by subtracting all soft and collinear limits. This means that for a phase-space 
point that lies very close to a singular point, the integrand is the square of the difference 
of two large numbers, namely the tree amplitude and its soft or collinear limit. In order to 
obtain the desired cancellation it is crucial to compute this difference in a numerically stable 
way, even if a certain invariant mass becomes very small. Thus, if the phase-space point 
is so close to a singular point that a straightforward evaluation of the amplitude becomes 
unstable, the amplitude is replaced by its (more stable) soft or collinear limit. We checked 
that the error introduced by this treatment is completely negligible. We also checked that our 
results are independent of the arbitrary parameters 5 and £ cut which have to be introduced 
in the subtraction method p3|. 


Another potential numerical problem is related to spurious singularities in the one-loop 
amplitudes. Besides the expected poles in the soft and collinear limits (which are avoided by 
the program since they lie in the three-jet region), the one-loop amplitudes have unphysical 
poles, i.e. poles with zero residue. Unfortunately, it is not possible to eliminate all these poles 
analytically, as long as the amplitude is expressed in terms of logarithms and dilogarithms 
multiplied by kinematic coefficients m ; this elimination is only possible if the amplitude is 
rewritten in terms of more general functions [^4]]. However, in the helicity formalism, one 


can simplify the (di)logarithmic coefficients to greatly alleviate the spurious poles We 


checked that the numerical evaluation of the matrix elements as given in refs. i 0 is stable, 
even for points that are quite close to a spurious pole, and that the probability for hitting an 
unstable point in the Monte Carlo integration is very small. Indeed, we had to evaluate close 
to a million points in a test run (corresponding to sub-percent statistical accuracy on the 
integral) in order to find one point that was ‘close’ to a particular spurious pole; at that one 
point the value of the vanishing denominator was still about an order of magnitude larger 
than where the numerical evaluation of the cross-section typically becomes unstable. 


3 Fixed-order Results 


We now present results for the four-jet fraction R A at next-to-leading order in a s . We use 
N c — 3 colors, Nf = 5 massless quarks, a strong coupling constant of a s (Mz ) = 0.118, a top 
mass of rri t = 175 GeV, a Z° mass and width of Mz = 91.187 GeV and Y z = 2.490 GeV, 
and a weak mixing angle of sin 2 9\y = 0.230 ||25|| . The numerical results given here are for 
a/s = Mz, but to the extent that contributions (II) and (III) can be neglected, R 4 depends 
essentially only on N c , Nf and a//?). We consider the E0, Durham |^6. 11| and Geneva (l^| 
jet algorithms. These cluster algorithms begin with a set of final-state particles (partons in 
the QCD calculation) and cluster the pair {i,j} with the smallest value of a dimensionless 
measure y l3 into a single “proto-jet”. The procedure is repeated until all the y l3 exceed the 
value of the jet resolution parameter y cut , at which point the proto-jets are declared to be 
jets. The algorithms differ in the measure y l3 used and/or in the rule used to assign a four- 
momentum pij to two clustered momenta Pi , pj. The same value of y cut in different schemes 
may sample quite different classes of events. For the reader’s convenience, we collect the 
definitions of the E0, Durham and Geneva schemes in Table []]. 

We start the presentation of the results with the E0 scheme. Fig. E§> shows the absolute 
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Tabic 1: Jet algorithm definitions 


Algorithm 

Uij 

Pij 

E0 

Durham 

Geneva 

(vi+PiY 

S 

0 min ( F 2, z?2\l-cos6y,- 

(^ + ^')(i, R5,) 

Pi + Pj 

Pi + Pj 

~ rmryn^ , n,j ) s 

8 F F 1 ^ cosdi i 

9 'A (Ei+EiY 2 


Table 2: EO algorithm 


Contribution to R 4 

Vent = 0.005 

Vent = 0.01 

Vcut = 0.03 

Born 

(2.60 ±0.02) • 10” 1 

(1.16 ±0.01) • 10” 1 

(1.79 ±0.01) • 10" 2 

a 

(2.43 ±0.08) • 10” 1 

(1.27 ±0.03) • 10' 1 

(2.42 ±0.05) • 10" 2 

b 

-(1.23 ±0.02) • 10' 1 

-(4.75 ± 0.04) • 10~ 2 

-(5.57 ±0.06) • 10" 3 

c 

-(4.06 ±0.02) • 10' 3 

-(1.83 ±0.01) • 10" 3 

-(2.93 ±0.01) • 10~ 4 

d+f 

—(1.13 ±0.18) • lO’ 2 

-(1.01 ±0.08) • lO' 2 

-(2.42 ±0.10) • 10~ 3 

e 

(1.42 ±0.01) • 10~ 2 

(5.45 ± 0.04) • 10’ 3 

(6.69 ±0.06) • 10~ 4 

II 

(1.66 ±0.28) • 10" 7 

(2.43 ±0.32) • 10” 7 

(1.88 ±0.18) • 10~ 7 

III 

-(1.18 ±0.01) • 10" 4 

-(7.53 ± 0.03) • 10' 5 

-(2.37 ±0.02) • 10” 5 

Full = R 4 

(3.79 ±0.08) • 10" 1 

(1.88 ±0.03) • 10" 1 

(3.46 ±0.05) • 10" 2 


value of the contributions of the different electroweak/color pieces to the four-jet fraction 
at y/s = Mz, as a function of y CVLt , setting the renormalization scale to n — Mz . Note 


that C[ b) + C 4 


(c) 


C 4 d> ± C 4 ' > and C 4 n) are negative. 


These curves are 


(from Table 

compared to preliminary SLD data points [H| which have been corrected for detector effects 
and hadronization. Obviously the comparison would benefit from a re-analysis using the 
full current Z° pole data samples. As expected, the subleading-color pieces are roughly 
10% of the corresponding leading-color contributions, reflecting the 1/N 2 suppression. This 
feature holds separately for the terms lacking and having an Nf factor. The contributions 
(11) and (III) are so small that we multiply them by a factor of 1000 and 10 respectively 
in the figure. Table presents the same results, for y cut E {0.005,0.01,0.03}, namely the 
coefficients (a s /27r) 3 0^/(1 ± 24 ) at y/s = Mz, including the statistical uncertainties from 
Monte Carlo integration. The ‘Born’ line gives the tree-level result (a s /2Ti) 2 B 4 /{1 ± —). 

Observable quantities calculated in QCD should be independent of the arbitrary renor¬ 
malization scale ye. However, the perturbative expansion is invariably truncated at a finite 
order, leading to a residual dependence of the result on /i. The tree-level /j dependence is 
much stronger for the four-jet rate than for the three-jet rate, because the former is propor¬ 
tional to a 2 instead of a s . The full ^-dependence of the NLO four-jet rate is given by 
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EO scheme 




/VVs 


Figure 2: (a) Absolute value of the contributions of the different electroweak/color pieces 

to the four-jet fraction at y/s = Mz for the EO scheme, i.e. (a s /2ir) 3 with 

x G {a, b, c, d, e, f, II, III}. We also show the Born and full one-loop prediction, and data from 
ref. [II]]. (b) Dependence of the tree-level (dashed line) and one-loop (solid line) prediction on the 
renormalization scale /r for y cut = 0.015. 
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Table 3: Durham algorithm 


Contribution to R 4 

2 /cut = 0.005 

Vent = 0.01 

Vent = 0.03 

Born 

(6.78 ±0.02) • 10~ 2 

(2.87 ±0.01) • 10~ 2 

(4.11 ±0.01) • 10" 3 

a 

(6.60 ±0.13) • 10- 2 

(3.03 ±0.06) • 10~ 2 

(4.23 ±0.07) • 10" 3 

b 

-(2.68 ±0.02) • 10" 2 

-(1.03 ±0.01) • 10" 2 

-(1.24 ±0.02) • 10" 3 

c 

-(1.27 ±0.01) • 10~ 3 

-(5.16 ±0.02) • 10~ 4 

-(6.94 ±0.02) • 10" 5 

d+f 

-(4.54 ±0.41) • 10~ 3 

-(2.50 ±0.09) • 10" 3 

-(3.67 ±0.45) • 10" 4 

e 

(2.93 ±0.02) • 10~ 3 

(1.14 ±0.01) • 10~ 3 

(1.43 ±0.01) • lO" 4 

11 

(2.28 ±0.20) • 10" 7 

(2.22 ±0.12) • lO” 7 

(9.06 ±0.39) • lO" 8 

111 

-(5.57 ±0.03) • lO' 5 

-(3.16 ±0.02) • 10" 5 

-(7.82 ±0.07) • 10~ 6 

Full = R 4 

(1.04 ±0.02) • 10' 1 

(4.70 ± 0.06) • lO’ 2 

(6.82 ±0.08) • 10” 3 


Table 4: Geneva algorithm 


Contribution to R 4 

2 /cut = 0.02 

2 /cut = 0.03 

2 /cut = 0.05 

Born 

(2.63 ±0.02) • 10 ” 1 

(1.50 ±0.01) • lO ’ 1 

(6.33 ±0.02) • 10 " 2 

a 

(1.16 ±0.05) • lO ' 1 

(8.91 ±0.25) • 10 ~ 2 

(4.90 ±0.14) • 10 " 2 

b 

-(1.37 ±0.02) • 10 ” 1 

-(6.99 ± 0.09) • 10 " 2 

-(2.51 ±0.03) • 10 ~ 2 

c 

-(7.78 ±0.12) • 10 " 3 

-(4.32 ±0.04) • 10 " 3 

-( 1.68 ± 0 . 02 ) • 10" 3 

d±f 

(6.90 ± 1.07) • 10 ' 3 

-( 1.10 ± 1 . 88 ) • 10” 3 

-(2.55 ±0.78) • 10 ~ 3 

e 

(1.44 ±0.02) • lO ' 2 

(7.58 ±0.08) • 10 ~ 3 

(2.83 ±0.03) • 10 " 3 

11 

(1.72 ±0.52) • 10 ” 7 

(2.89 ± 0.47) • lO ’ 7 

(2.53 ±0.35) • 10 ” 7 

111 

-(1.06 ±0.02) • lO ' 4 

-(7.86 ± 0.06) • 10 " 5 

-(4.91 ± 0.04) • 10 ” 5 

Full = R a 

(2.56 ±0.06) • 10 " 1 

(1.71 ±0.03) • lO ' 1 

(8.58 ±0.15) • 10 ” 2 


where a s (/x) is the two-loop running coupling, 


QL S (/i) 


W 


ots{M z ) 

w 



a s (M z ) (3i ln(w) \ 

7T Pq w ) 


7T 



( 12 ) 


with ft = i (fC A - IN,), ft = l(f C\ - (| C A + Cr)N,), C A = N c , C F = (JV= - 1)/(2JV C ). 
As expected, the strong renormalization-scale dependence of the tree-level result is reduced 
by the inclusion of the next-to-leading order contribution. Fig. ^b plots the /x-dependence 
of R 4 at tree-level and at one-loop for the EO scheme, at y cut = 0.015. 

The results for the Durham scheme are presented in Table [^, for the same values of y cut 
as in the EO scheme. Again, the subleading-color terms are of the expected size. 

The Geneva algorithm has the feature that the leading-order results, evaluated at /x = y/s, 
give a reasonable description of the data for large values of 2/cut, although the shape of 
the prediction is not quite correct, especially at small y cut ■ Also, the renormalization-scale 
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Geneva scheme 



Figure 3: NLO prediction for the four-jet rate using the Geneva algorithm for Nf = 5 and 
Nf = 8 . The theoretical bands have been obtained by varying the renormalization scale from 
< n < 2^fs and from < y < 3\/i. The data are from ref. [p~9|]. 


dependence is quite flat at moderate y cut . Finally, the dependence of the prediction on the 
number of light flavors Nf is reasonably large, at least in comparison with other algorithms 
(see Table [|). There is some interest in experimentally constraining Nf, in particular because 
a massless gluino would effectively shift the value of Nf by ANf = +3 in 0(ai%) four-jet 
distributions [13|, [T|. (At O(al) the effect is not simply given by ANf = +3, as is illustrated 
by the structure of the 0(a (?) results for the total e + e~ hadronic cross-section |27|.) Various 
authors have suggested that the existence of a light gluino is already in doubt [|L4|, |28| , 
Nevertheless, we would like to ask whether one can determine Nf with sufficient accuracy 
solely from the overall four-jet rate in the Geneva algorithm. In Fig. |3] we plot the NLO 
Geneva prediction as a function of y CVLt for Nf = 5 ( u , cl, s, c and b quarks) and Nf = 8 
( u,cl,s,c and b quarks, plus a massless gluino), where the bands represent the variation of 
y over the interval [^^/s,2y/s\ and 3^/s] respectively. These bands are compared to 

preliminary SLD data |TP[. The huge uncertainty for small values of y CVLt reflects the fact 
that the fixed-order prediction is not converging well for y cut < 0 . 02 , presumably due to 
large logarithms of l/y C ut■ This breakdown happens at larger y CVLt for Nf = 8 , since in this 
particular case Cf 1 is the dominant contribution to the one-loop correction, and it is further 
enhanced if Nf is increased from 5 to 8. 

As can be seen in Fig. |3] the data tend to favor Nf = 5, at least for 0.03 < y cnt < 0.04, 
however, the uncertainties coming from uncalculated higher-order terms are still too large 
to permit excluding light gluinos using this observable alone. 

Various angular distributions in four-jet events have been proposed to help separate the 
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relatively small contributions of four-quark final states from the dominant two-quark two- 
gluon final states |l2j . These distributions have been studied at leading order in a s in order 
to constrain iVj as well as the other color factors C A and C F m 0 0- The next-to- 
leading-order corrections to the distributions will be discussed elsewhere [TjJ, but they are 


remarkably small, given the size of the corrections to the overall four-jet rate. Unfortunately, 
in many cases the dependence on Nf is not that strong, such that a precise determination 
of Nf is difficult in the face of ha.dro ni za.tion uncertainties. 


4 Resummed Results 


The four-jet fraction declines rapidly at large y cut , and there is little data publicly available 
with which to compare our predictions for y cut > 0.07. On the other hand, at the kinematic 
limit y cut —>► 0 the QCD expansion parameter becomes a s L 2 , where L = ln(l/y cut ), and the 
NLO prediction would be improved if these large logarithms could be resummed. This is 
possible at leading order (LL) and next-to-leading order (NLL) in L in the Durham clustering 
algorithm because the phase space factorizes appropriately JTT]] . The NLL four-jet rate is 
then given by ||TT 


dNLL 

XL4 


= 2{A q (Q)} 2 \(r dqT q (Q,q)A 

. On 


Q 0 

[ dqT q (Q,q)A g (q) f dq (T g {q, q) A g (q') + T f (q') A f (q')) 
’ Q 0 J Q 0 


,(13) 


The NLL emission probabilities are 

r q (Q,q) = 
r g (Q,q) = 
r f(q) = 


2 C F a s (q) fQ _ 3 \ 
7T q \ q A) 
2C A a s {g ) f ln Q_ 11 N 


7 r q 

N f a s (q) 


q 12 


37T q 

and the Sudakov factors (probability of no emission) are 


(14) 


A q (Q) 
A g (Q) 
A f (Q) 


rQ 


= exp 


dqT q (Q,q)) , 
iQo J 

exp J Q dq\T q (Q,Q) + r /(?)] ) 1 

[A q (Q)Y 


A g (Q) 


(15) 


The Durham four-jet rate is an example of a quantity that can be resummed at leading 
and next-to-leading logarithmic order, but which does not exponentiate. The NLL results 
for such quantities do not include the proper renormalization-scale dependence of even the 


leading-log terms []30j]: Under a change of renormalization scale, a leading term a™L 2n varies 
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by ~ a™ +1 L 2n = a'” +1 Z/ 2< > n+1 )~ 2 , which is not contained in the NLL approximation. This is 
reflected in a relatively large ‘scale uncertainty’. Thus one should not rely on the resummed 
R 4 alone for a determination of . 

Indeed at finite values of y cut one should match the resummed results with the fixed- 
order results. For observables that exponentiate, a number of matching schemes have 
been defined [j3l| [|| — /^-matching, In R- matching, modified /i-matching and modified In R- 
matching. For i? 4 , the following matching scheme corresponds to R- matching: 


R 


R —match 


= i^ LL + 


a 


-V 


27 tJ 


(b 4 - 5 4 nll ) + (C 4 - C 4 nll ) 


1 + 


«< 


7T 


-1 


(16) 


where the ‘overlap’ terms I^ LL and C^ LL are defined by expanding i?^ LL out in powers of 
a s , in analogy to Eq. ( 0 )- A modified A-matching scheme could be defined by replacing 
L = ln(l/2/ cu t) by ln(2/hut — 2/mL + 1) i n -R 4 LL , where 7/ max is the maximum kinematic value of 
?/ cu t. This scheme would switch the resummed prediction over to the fixed-order prediction 
more quickly as y cut increases, and might therefore be more reliable at large y cut , but we 
have not yet implemented it. One could try to define an analog of In A-matching by 


7~>ln R —match 

4 


dNLL 

4 


b 4 

LL eXp 


a s 

2n 



C' 4 nll V 

5 4 nll JJ ’ 


(17) 


but Bf LL vanishes for y cnt ~ 0.01, so this approach fails. 

We evaluate the resummed _R^ LL using the two-loop formula (|l^) for the running coupling 
appearing in Eq. (|T4]). To evaluate the renormalization-scale dependence of i? 4 LL we make 
the substitution a s —> a s + /do ln(y 2 /s)a 2 /27T. In Fig. [| we show the resummed and matched 
prediction i?| J_match for the Durham algorithm, together with the tree-level and one-loop 
fixed-order predictions. In order to illustrate once more that the subleading-color terms are 
small we also show the leading-color one-loop result in Fig. f|. 

The agreement between theory and data is spectacularly good for the resummed and 
matched prediction. On the other hand, the ‘scale uncertainty’ in the prediction is still 
sizable. This is illustrated in Fig. [5] where the full one-loop and the resummed and matched 
results are shown as bands. These bands have been obtained by varying the renormalization 
scale from \Mz < y < 2Mz and \Mz < y < 3Mz respectively. (The large scale-dependence 
at large y cut in the resummed and matched prediction might be improved by a modified 
matching scheme.) 
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Durham scheme 



y cut 


Figure 4: The four-jet fraction for the Durham algorithm at y/s = Mz, illustrating the improve¬ 
ments to the Born term from adding successively the leading-color loop corrections, the subleading- 
color corrections, and the resummed corrections after matching. The data are from ref. [ 191. 
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Figure 5: Dependence on the renormalization scale of (a) the full one-loop prediction and (b) the 
resummed and matched result, for the Durham algorithm at \/s = Mz- 
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5 Conclusions 


In this article we presented the complete 0(a^) results for four-jet production in electron- 
positron annihilation. Generally, the NLO corrections are large and improve the agree¬ 
ment between theory and experiment considerably. The 1 /A^-suppressed correction terms 
are indeed smaller than the leading-color terms by the naive factor of ten or so. For the 
Durham algorithm, after the large logarithms of 1 /y cut have been resummed and the result 
is matched to the fixed-order prediction, and evaluated at the renormalization scale /i = y/s, 
theory agrees remarkably well with Z° data. Because the NLO corrections to the overall rate 
are so large, significant renormalization-scale dependence remains for both the fixed-order 
and resummed predictions, suggesting that there are still ~ 10 — 20% uncertainties from 
uncalculated higher-order corrections. More precise NLO predictions are possible for nor¬ 
malized four-jet distributions, for example the angles defined in ref. |12[, and will be reported 
elsewhere 


17 


Acknowledgement 

We thank Zvi Bern, Phil Burrows and David Kosower for valuable conversations and 
suggestions. 


Note added in proof. After we submitted this manuscript, Campbell, Glover and Miller re¬ 
ported on an independent calculation of the type (I) contributions to the one-loop virtual 
matrix elements for e + e _ —> qqgg |32|]. We have subsequently compared the virtual matrix 
elements used in this paper |S], |TU] to the results of refs. |32], , and we find agreement for 

both the four-quark and the two-quark-two-gluon final states. We thank J.M. Campbell and 
E.W.N. Glover for providing us with numerical results from refs. [32, |||. Also, Nagy and 


Trocsanyi |33|] have recently repeated our NLO calculation of li\ for the Durham, Geneva 
and E0 algorithms, using the matrix elements of refs. ||32| . ||. They obtain general agree¬ 
ment with our results, within statistical uncertainties. However, a relatively large difference 
(compared to the statistical errors) occurring at large y CVLt for the Geneva algorithm needs 
further investigation. 
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